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ABSTRACT 



o 
u 

o ' The mean color of globular clusters (GCs) in early- type galaxies is in general bluer 

than the integrated color of halo field stars in host galaxies. Metal-rich GCs often appear 
more associated with field stars than metal-poor GCs, yet show bluer colors than their 

jrt ■ host galaxy light. Motivated by the discovery of multiple stellar populations in Milky 

Way GCs, we present a new scenario in which the presence of second-generation (SG) 
stars in GCs is responsible for the color discrepancy between metal-rich GCs and field 
stars. The model assumes that the SG populations have enhanced helium abundance 
Q> ' as evidenced by observations, and it gives a good explanation of the bluer optical colors 

of metal-rich GCs than field stars as well as strong Balmer lines and blue UV colors 
of metal-rich GCs. Ours may be complementary to the recent scenario suggesting the 



f^ . difference in stellar mass functions (MFs) as an origin for the GC-to-star color offset. 



A quantitative comparison is given between the SG and MF models. 

Subject headings: globular clusters: general — stars: abundances — stars: evolution 
stars: horizontal-branch — stars: luminosity function, mass function 



1. INTRODUCTION 

Globular clusters (GCs) are the oldest stellar systems in galaxies and thus usually considered 
as the most representative tracer of early star formation in parent galaxies. Therefore, after the 
first discovery of bimodal color distribution of GCs in early-type galaxies (ETGs) together with the 
mean color offset (i.e., mean metallicity offset) between GC systems in ETGs and their host field 
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stars, the color discrepancy between GC systems and th eir host galaxies was generally interpreted 



as a different evolutionary h i story of these two systems (IForbes &: Fortd I2OOII: iBeaslev et al. 



20031: Forte et alJl2005l . l2007l : IPipino et alJl2007l : lLiu et alJl201ll ). Recently. iGoudfrooii k Kruiissen 



2002 



(|2013l ) reported that, even though the mean metallicity of red GCs, i.e., metal-rich GCs, and their 
field stars is almost identical, red GCs in ETGs show bluer colors than integrated colors of field 
stars in the host galaxies. In order to interpret this observation, they suggested that star clusters 
which initially formed with a "bottom-heavy" initial mass function (IMF) came to have a "top- 
heavy" IMF due to the long-term dynamical evolution, and their models reproduce the color offset 
between GCs and field stars of their host galaxies at the same metallicity (hereafter, in short, we 
will call this as the GC-to-star color offset). 

In this Letter, we suggest an alternative scenario that can also explain the GC-to-star color 
offset observed in ETGs. Our scenario is based on the presence of helium-enhanced stellar popula- 
tions that can reproduce extremely blue horizontal branch (HB) stars found in the Milky Way GCs 
(MWGCs) with multiple-generation stellar populations. As shown in IChung et al.l ( 201ll ). the hot 
HE stars from a helium-enhanced subpopulati on naturally reproduce the extremely blue integrated 



(FUV - V)o colors of metal-rich GCs in M87 (jSohn et al. 



2006 



Kaviraj et al.ll2007l ). and the effect 



of helium-enhanced subpopulations on the model (FUV — V)o color of metal-rich GCs is far greater 
than that of metal-poor GCs. We extend this result to the optical {g — z)o and {B — I)q colors 
and explain the GC-to-star color offset without applying IMF slope variations in the population 
models. 



2. POPULATION SYNTHESIS MODELS 



The models presented in this Letter are based on the Yonsei Evolutionary Population Synthesis 
(YEPS). The det ailed prescrip t ions f or the construction of the YEPS model are presented in the two 
recent papers of IChung et al.l (120111 ) and IChung et al.l (120131 ) . We have applied helium-enhanced 
subpopulations in the models to investigate the effect of second-generation (SG) populations with 
enhanced helium abundance on metal-rich stellar populations. 

Figure [T] illustrates the effect of helium-enhanced stellar populations on the morphology of 
HB stars with [Fe/H] = 0.0. We assume that the age of the helium-enhanced SG populations in 
the synthetic color magnitude diagrams (CMDs) is 1 Gyr younger than that of the first-generation 
(EG) population with Y = 0.23. In a single simulation, the number of HB stars in the EG and 
SG models is approximately 350 and 560, respectively. In addition, in order to avoid small number 
statistics in HB modeling, we have repeated the simulation 10 times to get the averaged number 
of HB stars in EG and SG populations. The larger number of HB stars in SG models is caused 
by their lower mean mass of HB stars that makes their lifetime longer. Since the mean mass 
of HB stars with normal helium abundance {Y = 0.23) is mainly controlled by metallicity and 
age, metal-rich and relatively young stellar populations show red HB types in the CMDs. These 



effects hold for the stellar population with Y = 0.28 at 11 Gyr. However, the effect of heliurrij 
on the morphology of HBs overtakes the effect of metallicity and age when the helium abundance 
reaches Y = 0.33. In the right panel of Figure [H although metallicity of the stellar population is 
as high as [Fe/H] = 0.0, the morphology of HBs of the stellar population with Y = 0.33 shows 
an extremely blue HE type. While the color difference of turn-off stars is < 0.2 mag between the 
Y = 0.23 and 0.33 populations, the mean color difference of HB stars is over ~ 2 mag in {g — z)q 
color. This means that the integrated optical colors of metal-rich stellar populations with helium- 
enhanced SG populations can be substantially bluer than that of models for the FG population 
with normal-helium abundances only. 

Using the SG population with enhanced helium abundance, we have constructed the model 
for GCs and field stars as follows. Firstly, we adopt Y = 0.33 for the helium abundance of the 
SG population and set the fraction of the SG populatio n for GCs to be 70 %. Th ese values are 



comparable to the simulated mean a bundance of heliuni (ILee et al.ll2005l : lJoo &: Ledl2013l ) and the 
observed fraction of SG populations (JD'Antona Sz Caloill2008l ) in the MWGCs with multiple stellar 
populations. One may argue that the assumed SG fraction in our model is rather too high, but 
note that the fraction of SG populations in some GCs like M 5, M13, M15, M53 , and NGC 6397 
is estimated to be more than 70 % based on the analysis of iD'Antona fc Caloil (J2008l ). Several 
lines of eviden ce suggest that 7 - 8 % of stars in the inner halo were originated from SG stars in 



massive GCs (IMartell et al. 
with Y 



2011 



Therefore, we have also adopted 10 % of the SG population 
0.33 for field stars in host galaxies because the MW halo might be regarded as the local 
counterpart of ETGs. S econd ly, we assume the age gap between two generations as 1 Gyia. As 
discussed in I Joo &: Led (120131 ) . the age difference between the FG and SG stars for the MWGCs 
with an enhanced helium abundance is between 0.3 — 1.7 Gyr. For all models with SG populations, 
we have adopted the standard Salpeter IMF, ^(m) oc m~", with a slope of a = 2.35. In order to 
examine the effect of the IMF slope change on the integrated colors of simple stellar populations, 
we have also simulated models with different IMF slopes of a = 1.7 and 3.0 for FG populations. 
We note that our simulation is designed for both GCs and field stars in ETGs, and thus applying 
the model to other types of galaxies should be done with caution. Table [1] summarizes the input 
parameters of models used in this Letter. 



^We consider the helium content as the third parameter of determining HB types in MWGCs with multiple 
generation stellar populations, and the difference in helium conte nt among MWGCs is attributed to the d ifferent 
formation efficiency of SG stars with enhanced helium abundance I D'Antona fc Caloil l2008l : I Joo fc Leell2013l l. 



The significant fraction of ETGs shows recent star formation (RSF) signatures (|Yi et al 



2011 



However, the 
most ETGs with strong bimodal GC color distributions show UV-upturn phenomenon with no signatures of RSF, 



and the effect of RSF on the integrated colors of ETGs is almost negligible in optical colors (see Fig 4 of lYi et al 
201 ll ). Therefore, we did not include RSF in our simulation for host galaxies. 



3. RESULTS 

Upper panels of Figure [2] show the effect of the SG populations on the integrated optical 
colors of {g — z)q and {B — I)q, as well as (FUV — V)o. For the comparison of models with 
the observation, we overplot ob served colors of red GC s and field stars in the left and middle 
panels using data in Table 1 of iGoudfrooij &: Kruijssenl (120131 ) . In our simulations, we assume 
the metallicity of metal-rich GCs as [Fe/H] = O.Qj becaus e spectroscopy of typical metal-rich GCs 
suggests the solar metallicity in various observations (e.g.. lCohen et al.ll2003l ). We also presume the 
observed metallicity of field stars to be [Fe/H] =0.0. In the metal-rich regime, blue and extremely 
blue HB stars from helium-enhanced populations make the integrated color of stellar populations 
bluer. The strongest effect of hot HB stars can be observed in the (FUV — V)q color (the right 
panel in the top row). This effect is not limited to (FUV — V)o and makes {g — z)q and {B — I)q 
of the SG population, on average, 0.24 and 0.35 mag bluer than those of the FG population (see 
the gap between blue and red dotted lines on the gray shades in the top panels). The integrated 
colors of {g — z)o and {B — I)q for the metal-rich GC model with the 70 % of SG populations show, 
respectively, 0.14 mag and 0.21 mag bluer colors than the model for field stars of the host galaxy 
with the 10 % of SG populations at [Fe/H] = 0.0. Considering the r eported GC-to-star color offset s 
of ETGs ranging from 0.14 to 0.23 for {g-z)o and 0.21 for {B - I)o dCoudfrooii fc KruiisseJ2013l ). 
our simulation for [Fe/H] = 0.0 well explains GC-to-star color offsets observed in the ETGs without 
applying the IMF variation. 

We have also tested the effect of the IMF variation on the optical colors in the bottom panels of 
Figure[2j We assume two different models with the bottom-heavy IMF (a = 3.0) and the top-heavy 
IMF (a = 1.7) compared to the standard Salpeter IMF (a = 2.35). The color offsets caused by 
the variation in the IMF slope from bottom to top heavy, at the same metallicity ([Fe/H] = 0.0), 
are 0.07 mag and 0.06 mag for {g — z)q and {B — I)q, respectively. These offsets are not enough 
to account for the observed GC-to-star offsets in (g — z)o and {B — I)q. Moreover, the model 
with variations in the IMF hardly reproduces the observed (FUV — V)q colors of GCs in M87 (see 
the right most panels). It is also no teworthy that the present-day mass function ( MF) of clusters 
deviates significantly from the IMF (JBaumgardt &: Makindl2003l : lPaust et al.ll2010l ). However, the 
effect of the present-day MF on the integrated color of GCs is very small. The color offsets caused 
by the present-day MF (a = 1.7) are less than 0.01 mag in both {g — z)q and {B — I)q when the 
[Fe/H] is 0.0. This result implies that, in order to explain GC-to-star color offsets observed in 
various colors simultaneously, the inclusion of a SG population in the model is more likely than 
adopting different IMF slopes in the model. 

Since the absorption strength of H/3 is very sensitive to the hot temperature (Teff ~ 10, 000 K) 
of blue and extremely blue HB stars, it is important to check the effect of the SG population on H/3. 



Note that if the metaUicity of red GCs decreases, the effect of SG populations on the optical colors significantly 
reduces. Therefore, our simulation presented here can be better applied to metal-rich GCs whose metallicity is greater 
than [Fe/H] = -0.2. 



Figure [3] compares our model with observations of NGC 4472 and NGC 1407. The upper left panel 
presents the models for FG populations with ages from 1 to 14 Gyr and the model for the 11-Gyr 
SG population. Similar to the integrated colors, when [Fe/H] is 0.0, H/3 of the model with the 
SG population for 11 Gyr shows 1.05 A stronger absorption than the model of the FG population 
for 12 Gyr. This absorption strength is equivalent to the 3 Gyr model of the FG populations. As 
shown in the right top panel, the change of the IMF slope from top to bottom-heavy decreases the 
strength of H/3 for the metal-rich population of [Fe/H] = 0.0 by 0.11 A. This is only a 10 % change 
with respect to the effect of the SG populations. 

In the bottom panel s of Figure [3l we compare ou r models with the GCs in NGC 4472 and 
NGC 1407 (data are from lGoudfrooii &: Kruiissenll2013l ). The arithmetic mean of observed II/3 and 
[MgFe]' are also presented. Given the observational errors, both models for the FG and the SG 
population bracket the observed GCs in NGC 4472 and NGC 1407. The overall GC absorption 
indices of NGC 4472 and NGC 1407 are better explained by our model with SG populations because 
the inclusion of the SG populations in the model slightly enhances the strength of II/3 by 0.4 A 
when the mean [Fe/H] of the metal-rich GCs is 0.0. The H/3 difference between two observed 
quantities of NGC 4472 differs at ~ 1.96a confidence level. Although the confidence level of the 
H/3 difference in NGC 1407 is lower (~ 1.02o-), this is not inconsistent with the result of NGC 4472. 
In general, the age-dating of GCs in ETGs based on the model only with FG populations gives 
systematically younger ages for metal-rich GCs because of the sca ttered distributi on of metal-rich 



GCs. For example, the age of GCs in NGC 4472 (jCohen et al.ll2003l ^ and NGC 5128 (jWoodlev et al. 



20 id ) become younger with increasing metallicity. However, the inclusion of a SG population in the 
model reduces the age-metallicity gradient and reproduce simultaneously integrated optical colors 
and absorption indices of GCs in ETGs without any further fine tuning. This result reinforces 
that the helium-enhanced SG population is a very plausible candidate for the subpopulation in the 
metal-rich GCs in ETGs. 



4. DISCUSSION 



We have reproduced the observed GC-to-star color offsets found in massive ETGs using helium- 
enhanced SG stellar populations. Our models presented here suggest that some 70 % of helium- 
enhanced SG population explains GC-to-star color offsets observed in the ETGs without applying 
the variation in the IMF slope. In addition, the general trend between H/3 and metallicity found 
in metal-rich GCs can also be explained by our models with S G populations. T hese results are 



consistent with our previous models for metal-rich GCs in M87 (jChung et al.ll201ll ). which showed 



that the inclusion of helium-enhanced subpopulations in the models can provide the most reasonable 
explanation for the extremely blue (FUV — V)q color. 

Our result may shed some light on the long-standing problem of the disagreement between GCs 
and filed stars in their metallicities. The discrepancy in metallicities are twofold: (a) the metallicity 
distribution function (MDF) morphologies and (b) the mean metallicity values. lYoon et al.l ()201ll ) 



delved into the difference in MDF morphology and, based on their theoretical color-to-metallicity 
conversion relations, suggested that inferred GC MDFs are remarkably similar to the MDFs of 
resolved field stars in nearby ETGs. As for the mean metallicity values of GCs and stars, however, 
the discrep a ncy r emains in the sense that GCs are bluer than field stars within the same galaxy. 
Yoon et al.1 (|201ll ) attributed the remaining GC-to-star metallicity offset to the difference in their 
details of chemical evolution processes, proposing a scenario in which GCs are the relics of vigorous 
star formation early on, and thus trace selectively the major mode of star formation. That is, GC 
formation was less prolonged than field-star formation. As a consequence, the chemical enrichment 
of GC systems has ceased earlier than that of the field stars, leaving GCs that are metal-poorer than 
field stars. Our SG scenario predicts that metallicity of metal-rich GCs is the same as that of field 
stars, further reducing the degree of the m etallicity difference between GCs and stars. Our results, 
combined with those of lYoon et al.l (|201ll ). indicate that GC systems and their parent galaxies have 
shared a more common evolutionary history than previously believed. 



Since the investigation by Ivan Dokkum fc ConroyI (|201Cl ) , the idea of a bottom-heavy IMF for 



the massive ETGs has received much attention in the curren t literature (e.g-JCaPDellari et al. 



2012 



van Dokkum &: Conrovl 120121 : IConrov &: van DokkumI l2012l : IZaritskv et al.l |2013| ) . Their interpre- 



tations are based on the strengths of the Na I doublet at A 8183, 8195 A and the Wing-Ford band 
at A 9916 A that become stronger when low mass stars are abundant. However, it is well known 
that the strength of spectroscopic absorption line is controlled not only by the gravity but also by 
the temperature and chemical composition. Therefore, if field stars in ETGs are contaminated by 
the SG populations from massive GCs, the strong Na I doublet of ETGs may also be caused by SG 
populations that generally contain enhanced N a (jGratton et al.l 12004 : ICarretta et al.l |2005J, l2009l ) . 
The Figure 2b of Ivan Dokkum &: ConroyI (120101 ) gives some hints for the effect of the enhanced Na 
element because their result based on Na I doublet prefers steeper IMF slope (bottom-heavy IMF) 
than that based on the Wing-Ford b and measurement . In addition, the Salpet er-like IMF of ETGs 
are reported based on spectroscopy (jSmith et al.ll2012l ) and stellar dynamics (JThomas et al.ll201ll ) 
of ETGs. We will investigate this issue in a separate paper. 
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Fig. 1. — Effect of helium abundance on the morphology of HB stars in the color- magnitude 
diagram of simple stellar populations. From left to right panels, the helium abundance of stellar 
populations in each panel is 1" = 0.23 (red), 0.28 (green), and 0.33 (blue), respectively. The age of 
the helium-enhanced SG populations (y = 0.28 and 0.33) is 1 Gyr younger than that of the FG 
population with normal helium abundance {Y = 0.23). For the comparison with normal-helium 
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Fig. 3. — {top row) Effects of the SG population and the IMF slope on the integrated absorp- 
tion indices of the UP and the [MgFe]' (= ^/Mgb x (0.72Fe5270 + 0.28Fe5335) A). All black hnes 
indicate 100 % FG models. Red and blue lines in the left panel highlight 100 % FG and 100% 
SG models at the age of 12 and 11 Gyr, respectively. The red and blue lines in the right panel 
indicate different IMF slopes with the a = 1.7 and 3.0, respectively. The age and the metallicity 
are indicated in the plot, (bottom row) Comparison of the models presented in the Figure [2] with 
GCs in NGC 4472 and NGC 1407. Blue and red open circles are blue and red GCs in NGC 4472 
and NGC 1407. Typical error bars for each GC system are indicated. Filled blue and red squares 
indicate arithmetic mean of blue and red GCs in each galaxy NGC 4472 and NGC 1407. Black 
triangles are the integrated absorption index of stars in the host galaxies. 
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Table 1. INPUT PARAMETERS ADOPTED IN THE SECOND-GENERATION AND 

COMPARISON MODELS 



Parameters 



Models for GCs 



Models for host galaxy field stars Comparison models 



Slope of Salpeter initial mass function, a 2.35 

a-elements enhancement, [a/Fe] 0.3 

HB mass dispersion, ctm (^0) 0.015 

Reimers' mass-loss parameter, rj 0.63 

Absolute age, t (Gyr) 11.0 and 12.0 

Helium enhanced population ratio, Y (%) 0.23 (30%) + 0.33 (70%) 



2.35 


1.7 and 3.0 


0.3 


0.3 


0.015 


0.015 


0.63 


0.63 


11.0 and 12.0 


12.0 


90%) + 0.33 (10%) 


0.23 (100%) 



